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ABSTRACT 

The very high energy (VHE) 7-ray spectra of some TeV Blazars, after being 
corrected for absorption in the extragalactic background hght (EBL), appear un- 
usually hard, which poses challenges to conventional acceleration and emission 
models. We investigate the parameter space that allows the production of such 
hard TeV spectra within time-dependent leptonic models, both for synchrotron 
self-Compton (SSC) and external Compton (EC) scenarios. In the context of 
interpretation of very hard 7-ray spectra, time-dependent considerations become 
crucial because even extremely hard, initial electron distributions can be signif- 
icantly deformed due to radiative energy losses. We show that very steep VHE 
spectra can be avoided if adiabatic losses are taken into account. Another way to 
keep extremely hard electron distributions in the presence of radiative losses, is 
to assume stochastic acceleration models that naturally lead to steady-state rel- 
ativistic, Maxwellian-type particle distributions. We demonstrate that in either 
case leptonic models can reproduce TeV spectra as hard as dN/dE^ oc E^. 
Unfortunately this limits, to a large extend, the potential of extracting EBL from 
7-ray observations of blazars. 
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Subject headings: BL Lacertae objects: general - diffuse radiation - gamma-rays: 
observations - gamma-rays: theory. 



Introduction 



Blazars constitute a sub-class of Active Galactic Nuclei (AGN) characterized by broad- 
band (from radio to VHE 7-rays), non-ther mal emission produced in relativistic jets pointing 
close to the line of sight to the observer (jUrry &: Padovanilll995[ ). The highly variable lu- 
minosity of Blazars (which often exhibits 2 peaks) are commonly interpreted in terms of a 
synchrotron-inverse Compton origin. 

In synchrotron self-Compton (SSC) models, the X-ray emission is usually attributed 
to the synchrotron radiation of relativistic electrons. The Compton up-scattering of these 
synchrotron photo ns by the same electron population then prod uces the high energy 7- 



ray radiation (e.g., iMaraschi et al.lll992l : iBloom &: Marcher! 119961 ) . Under specific circum- 



stances, the target radiation field for inverse Compton upscattering can be dominated by 



external photons, leadin g to so-called external Compton (EC) models (e.g., iDermer et al. 
I993I : ISikora et al.lll994j ). In general, these leptonic models have been relatively successful in 
describing the observed SED of Blazars. 

The recent detections of VHE 7-rays from Blazars with redshift z > 0.1 (in particular, 
lES 1101-232 at ^ = 0.186 and lES 0229+200 at z = 0.139), however, poses challenges to 
the conventional leptonic interpretation. VHE 7-rays emitted by such distant objects arrive 
after significant absorption caused by thei r interactions with extrag alactic background light 



(EBL) via the process 77 — )■ e+e (e.g., iGould fc Schredeii 119671 ). Reconstruction of the 



absorption-corrected intrinsic VHE 7-ray spectra based on state-of-the-art EBL models then 
yields unusually hard VHE source spectra, that are difficult to account for with standard 
inverse Compton assumption. 

One characteristic case concerns the distant (at z = 0.186) Blazar l ES 1101-232, de- 



tected at VH E 7-ray energies by the H.E.S.S. array of Cherenkov telescopes ( lAharonian et al. 



2006|, |2007aJ) . When corrected for absorption by the EBL, the VHE 7-ray data result in very 
hard intrinsic spectra, with a peak in the SED above 3 TeV and a photon index F < 1.5. 
A similar behavior has also been detected in the TeV Blazar lES 0229+200 (at z=0.139) 
( lAharonian et al.ll2007bl ). Though there is a non-negligible uncertainty in the EBL flux, the 
i ntrinsic spectra are unus ually hard even when one considers the lowest levels of the EBL 
( iFranceschini et al.l 120081 ). Other models predicting higher EBL flux would lead to even 
harder (intrinsic) photon indices close to 1 (e.g., Stecker & Scully 2008). We note that a 



- 3 - 



recent analysis of Fermi LAT data for the nearby TeV Bl azar Mkn 501 indicat es a hard 7-ray 
spectrum (F close to 1) at lower (10-200 GeV) energies (INeronov et al.ll201ll ). If confirmed, 
this would be strong evidence for unusually hard 7-ray spectra independent of questions 
related to the level of EBL. On the other hand and apart from the challenges arising for 
inverse Compton interpretations, the observed hard VHE spectra obviously carry important 
information about the level of the EBL, and thus a deep understanding of the mechanisms 
acting within these sources becomes now even more critical. 

The "simplest" way to overcome the problem is to assume t hat there is n o absorption. 
In fact, this is possible in Lorentz invariance violation scenarios (IKifund 119991) . We shoul d 
note however that this effect is likely to be true only above 2 TeV (IStecker fc Glashowll200ll ). 
whereas the hard spectra problem we face in the case of distant Blazars is relevant to sub- 
TeV energies as well. Anot her non-standard mech anism to avoid severe absorption in the 
EBL has been suggested by iDe Angelis et al.l l2009l who proposed that 7-ray photons could 
oscillate into a new very light axion-like particle close enough to the source and be converted 
ba ck before reaching the Earth. To some extend a similar idea was recently suggested 
by lEssey et al.l 1201 ll who suggested that the 7-rays from Blazars may be dominated by 
secondary 7-rays produced along the line of sight by the interactions of cosmic rays protons 
with background photons. While the first scenario would require the existence of exotic 
particles, the second needs extraordinary low magnetic fields of the order of lO^^^G. 

In more standard-type astrophysical scenarios, the formation of hard 7-ray spectra is 
related to the production and absorption processes. Photon-photon absorption could in 
principle result in arbitrarily hard spectra provided that the 7-rays pass through a hot 
photon gas with a narrow distribution such that E^eo » nieC^. In this case, due to the 
reduction of the cross-section, the source becomes optically thick at lower energies and thin 
at higher energies, thus leading to formation of hard intrinsic spectra. ( lAharonian et al. 
20081 : IZacharopoulou et"aDl201l[ l 



If one relates the hard 7-ray spectra to the production process, then this implies corre- 
spondingly hard parent particle distributions. Outside standard leptonic models, a number 
of alternative explanat ions have been explored in the literature. In analogy to pulsar winds, 
haronian et al.l ( 120021 ) for example have analyzed the implications of a cold ultra-relativistic 
outflow that initially (close to the black hole) very high speeds. In such a case, 

up-scattering of ambient photons can yield sharp pile-up features in the intrinsic source spec- 
tra. However, very high bulk Lorentz factors would be needed (F;, ~ 10^) and it seems not 
clear whether such a scenario can be applied to Blazars. On the other hand, if Blazar jets 
would remain highly relativistic out to kpc-scales (Ff, ~ 10) and able to accelerate particles, 
a hard (slowly variable) VHE emission component could perhaps be produced by Compton 
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up-scattering of CMB photons (iBottcher et al.ll2008l ). 



In order to produce hard 7-ray spectra within standard leptonic, synchrotron-inverse 
Compton scenarios, hard electron energy distributions are required. Ahhough standard 
shock acceleration theories, both in the non-relativistic and relativistic regime, predict quite 
broad, n{E) oc -E'~^-type, electron energy distrib utions, there are non - conventional realiza - 
tions which could give rise to very hard spectra (IDerishev et al.l l2003l : IStecker et al.l 120071 ) . 
On a more phenomenological level, iKatarzyhski et al.l ( 120071 ) have shown that the presence 
of an energetic power-law electron distribution with a high value of the minimum cut-off 
energy can lead to a hard TeV spectrum. In general, however, injection of a hard electron 
distribution is not a sufficient condition as electrons are expected to quickly lose their energy 
due to radiative cooling and thereby develop a standard n{E) cx E'"^ form below the initial 
cut-off energy. In order to avoid synchrotron cooling , one thus needs to a ssume unusually 
small values for the magnetic field within the source (ITavecchio et al.ll2009l ). 



To some extent, this situation can be avoided if one invokes adiabatic losses. This is 
demonstrated below by means of a time- dependent investigation assuming dominant adia- 
batic energy losses. As a second alternative, we discuss pile-up (Maxwellian-type) electron 
distributions, that can be formed in stochastic acceleration scenarios. As these distributions 
are steady-state solutions with radiative (synchrotron or Thomson) losses already included, 
there is no need to avoid these losses. Maxwellian-type electron distributions provide an 
interesting explanation for the very hard TeV components as their radiation spectra share 
many characteristics with the (hardest possible) mono-energetic distributions. 

It is obviously important to explore the strengths and limitations of such explanations 
in more details, both theoretically and observationally, in order to understand whether there 
is a need to invoke more exotic scenarios. 

In the present work we explore the conditions under which a narrow, energetic particle 
distribution is able to successfully account for the hard VHE source spectra. To this end, we 
examine different electron distributions within the context of standard leptonic models, i.e. 
the one-zone SSC and the external Compton scenario. The paper is structured as following: 
The requirements for quasi-stationary SSC solutions are analyzed in Sect. 2. Apart from 
narrow, power-law-type electron distributions, quasi-Maxwellian distributions are examined. 
A time-dependent generalization including adiabatic losses is explored in Sect. 3. Section 4 
discusses the possibilities within an external Compton approach. 
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2. Stationary SSC with an energetic electron distribution 

Within a stationary SSC approach, the hardest possible (extended) VHE spectrum 
is approximately F^, cx u^^^, where Fi, = dF/du is the spectral flux (differential flux per 
frequency band). This has a simple explanation: The emitted synchrotron spectrum of a 
single electron with Lorentz factor 7 in a magnetic fleld B, averaged over the particle's 
orbit, obeys oc G{x), where G{x) is a dimensionless function with x = u/uc and Uc = 

87^6-8 sin a/ (47rmec). For x ^ 1, the functional dependence o f G{x) is well approximate d 



by G{x) oc x^/'^, while for x ^ 1 one has G{x) oc x^^^e"^ (e.g., iRybicki fc Lightmanlll979l ). 
Hence, at low frequencies v <^ z/^, the synchrotron spectrum follows oc z/^/^. Compton 
up-scattering of such a photon spectrum in the Thomson regime by a very energetic, narrow 
electron distribution will preserve this dependence and therefore yield a VHE spectral wing 
as hard as F^, oc u^^^ (see below). 



2.1. Power-law electron distribution with high low-energy cut-off 

A homogeneous SSC scenario with a high value for the low- energy cut-off of the non 



thermal electron distribution has consequently been proposed by iKatarzyhski et al.l (120071 ) 
to overcome the problem of the Klein-Nishina (KN) suppression of the cross-section at high 
energies and to reproduce VHE spectra as hard as 1/3. Let us assume that the electron 
population follows a power-law distribution of index p between the low- and high-energy 
cut-offs 

as often used in modeling the Blazar spectra. Here, prime quantities refer to the blob rest 
frame and unprimed to the observer's frame. Taking relativistic Doppler boosting {6) into 
account, the observed synchrotron flux from an optically thin source at distance d^ is given 
by the integral of N'^{'~f')d'~f' tim es the single particle emissivity j'(z/',7') over the volume 



element and all energies 7' (e.g., iBegelman et al.lll984j ). i.e 



^r = |-/ I fWn')K{i)didv'. (2) 



The above expression yields the common power-law of index a = (p — l)/2 between the 
frequency limits i/min oc 5(57^;^) and z/max oc (5(i?7^ax)- Below and above those limits, the 
electrons with energy around the minimum and maximum cut-off dominate, and thus the 
spectrum approximately exhibits a slope Fy oc v^l"^ for v < z^min, and an exponential cut-off 
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for v > iy„ 



I.e. 



oc 



p-1 

U 2 



(3) 



l/2„-i. 



The hard 1/3-slope appears in the VHE range of EC 7-rays when the synchrotron 
photons are up-scattered to higher energies by the electron population given by equation ([1]) 
with a high 7min and provided that the Thomson regime applies. Obviously, it will be 
significantly softer in the KN regime. In any case, however, there exists a characteristic 
energy below which the Compton spectrum mimics the behavior of the synchrotron spectrum 
oc z/i/3. 

Note that the inverse Compton-scattered spectrum of a monochromatic photon field by 
mo no-energetic electrons appro ximately follows, at low up-scattered photon energies, F^ oc v 
(cf. iBlumenthal fc Gouldlll970l ). Thus, any photon field which is softer (flatter) than F^^ oc v 
will dominate the lower-energy part of the up-scattered emission and thus, in the standard 
SSC scenario the 1/3- VHE slope (the 4/3-slope in the vF^ representation) is the hardest 
that can be achieved. 

An exemption to this may occur if the magnetic field in the source would be fully tur- 



^^^^ . 

{ Medvedev 


2006; 


Derishev 


2007: 



2OIOI ). which will then be reflected to low-energy part of the Compton component. 



The "critical Compton energy" is usually emin — (^Tmml^Tmm)) where h = {B / Bcr)meC^ , 
Bcr = mlc^/{eh), except when the deep KN regime applies, i.e., when up-scattering of the 
minimum synchrotron photons by the minimum energy electrons occurs in the KN regime so 
that |&7min > 1- If the latter applies, then the corresponding energy below which one can see 



the hard 1/3-slope is, as expected, 7n 



and it approximately corresponds to the peak 



of the emitted luminosity for any powe r-law electron index (s ee Fig. [T]). In the KN regime, 
the peak appears especially sharp (e.g., iTavecchio et al.lll998l ). and the Compton flux has a 
strong inverse dependence on the value of 7min- For example, for the realization presented 
in Fig. (11]), the emissivity in this regime roughly scales as oc 7~fn^, so that slight changes 
in 7min can lead to significant variations in the amplitude of the Compton peak flux. On 
the other hand, as long as p < 3 (positive synchrotron slope in a uF^, representation) the 
synchrotron peak luminosity would remain approximately constant. 



A power-law electron distribution with a high low-energy cutoff has been used in ITavecchio et al 



( I2OO9I ) in order to reproduce the SED of the blazar lES 0229-1-200 within a stationary SSC 
approach. The high value of 7min ~ 10^ then ensures the hard Compton part of the spec- 
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trum with 1/3-slope is in the TeV range. The generic difficulty for such an approach is 
that an energetic electron distribution is expected to quickly develop a 7~^-tail below 7min 
due to synchrotron cooling, thereby making the C ompton VHE spectrum softer (see Fig. H]). 
To overcome this problem, iTavecchio et al.l (120091 ) suggested an unusually low value for the 
magnetic field, B ~ (10~^ — 10~^) G, that would allow the electron distribution to remain 
essentially unchanged on timescales of up to a few years. Obviously, one would then not 
expect to observe significant variability on shorter timescales. We note however, that this 
requirement could be relaxed if one assumes that the detected 7-ray signal is a superposition 
of short flares which can not be detected individually. Arguments based on magnetic flux 
conservation naively suggest that the magnetic field value, when scaled from the black hole 
region to the emission site, should be at least one or two orders of magnitude higher so that 
one would need to destroy magnetic flux for such a scenario to work. On the other hand, 
a narrow but very energetic electron distribution in combination with such low magnetic 
field strengths implies a strong deviation from equipartition, thereby obviously facilitating 
an expansion of the source. 



2.2. Relativistic Maxwellian electron distribution 



As far as a narrow energetic particle distribution is concerned, a relativistic Maxwellian 
may come as a more natural representation. Such an electron distribution can be the outcome 
of a stochastic acceleration process (e.g., 2nd order Fermi) that is balanced by synchrotron 
(and/or Compton) energy losses, or in general any ene rgy loss mechanis n i that exhibits a 



quadratic dependence on th e particle energy (see e.g., ISchlickeiseii Il985l : lAharonian et al. 
19861 : iHenri fc Pelletiei]ll99l[ Istawarz fc Petrosian|[2008h . 



Consider for illustration the Fokker-Planck diffusion equation which describes the sta- 
tionary distribution function f{p) of electrons that are being accelerated by, e.g., scattering 
off randomly moving Alfven waves in an isotropic turbulent medium. 



l_d_ 

p2 Qp 



+ 



1 d 



df{p) \ 
dp ) ' p^ dp 



(4) 



where Dp is the momentum-space diffusion coefficient. Particle escape is neglected in eq. (jl]), 
as the timescale for synchrotron cooling is expected to be much smaller than the one for 
electron escape. 



For scattering off Alfven waves, one has Dp 



wit h Va 



the A lfven speed and t = X/c oc ap > 0, the mean scattering time (e.g., iRieger et al. 
20071 ). If the turbulent wave spectrum W{k) cx k~'^ is assumed to be Kolmogorov-type 
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(g ~ 5/3) or Krai chnan- type (g = 3/2), the momentum-dependence becomes = 1/3 and 
Op = 1/2, respectively. Bohm-type diffusion, on the other hand, would imply Op = 1, while 
hard-sphere scattering is described by ap = 0. Note, however, that if one considers electron 
acceleration by resonan t Langmuir waves, even Dp = const {ap = 2) may become possible 



(lAharonian et al.lll986l ). 



The synchrotron energy losses that appear in the second term of Eq. (jlj) are 

f = = -^K/m^c^)f/V (5) 

In the 7-parameter space, the solution of Eq. (jl]) becomes a relativistic Maxwell-like 
function 

/(7) = Ve-(-)'''"' , (6) 

{ttp 7^—1) with 

7c =( ^ 1 [mec) , (7) 

and constant A to be defined by the initial conditions. Note that this is a steady-state 
solution already including radiative losses and there is no need to invoke extreme values for 
the magnetic field. The critical Lorentz factor 7c approximately corresponds to the energy 
at which acceleration on timescale 

\ 2 

c 



is balanced by (synchrotron) cooling on timescale tcooi = l/[/^sP]- Depending on the choice 
of parameters, a relatively large range of values for 7c is possible and thus, cut-off energies 
of the order of 7c ~ 10^ may well be achieved. Consider, for example, Bohm-type diffusion 
with r = rjrg/c, Vg = •yrrieC^ / (eB) the electron gyro-radius and r] > 1. Using tacc = ^cooh the 
maximum electron Lorentz factor becomes 7c — 10^ (f^i/O.Olc) (1 Q / BY^'^r]^^^'^ . 

The synchrotron spectrum that arises from a Maxwell-like electron distribution is dom- 
inated by the emission of electrons with 7c (Fig. [3l). It exhibits the characteristic 1/3-slope 
up to the corresponding "synchrotron cut-off frequency" hv^^^ ~ (567c where h = B/B^r and 
Bcr = m?c^/eh. Thus the Compton spectrum is very similar to the one resulting from a 
narrow power-law if one chooses a value for the cut-off energy close to the minimum elec- 
tron energy of the power-law distribution. The peak of the Compton flux then contains 
information for the cut-off energy as t'pgak 7c- 

Note that for an electron distribution of the form of eq. ([6]) that exhibits an exponential 
cutoff oc exp[— (7/7c)^], the corresponding cut-off in the synchrotron spectrum appears much 
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smoother, oc exp[— (zz/z/*)''/'^'^"'"^)] ( Fritzlll989l : IZirakashvili fc AharonianI 120071 ). The position 
of the synchrotron peak flux, z/p, is then also dependent on /3, and one can show that for /3 = 1 
(or Op = in the previous notation) an important factor ~ 10 arises, so that z/p = d.buc, 
whereas for /3 = 3 the synchrotron peak corresponds approximately to the electron cut-off 
as z/p = 1.2z/c (e.g.. Fig. [3]). 



3. Time-dependent case - expansion of the source 

Expansion of the source could change the conclusions drawn above. In particular, if one 
assumes a very low magnetic field such that synchrotron losses are negligible, then adiabatic 
losses may become important and alter the electron distribution. In this section, we examine 
the behavior of the system for a power-law electron distribution with a high value of the low- 
energy cut-off discussed above. For simphcity, we consider a spherical source that expands 
with a constant velocity u, 

R{t) =Ro + u{t - to) . (9) 
The relativistic electron population will be affected by synchrotron losses, 

at DTrmgC 
and by adiabatic losses (e.g., Longair 1982), 

rf7 kit) u 

As the emission region expands, the magnetic field decreases. We consider a scaling B cx 
(l/R)"^ oc (1/t)™ with 1 < m < 2 to study the evolution of the system. The limiting 
value m = 2 corresponds to conservation of magnetic flux for the longitudinal component, 
whereas m = 1 holds for the perpendicular component. (Note that for m = 1 the ratio of 
the electrons' energy density to the magnetic field energy density remains constant). Which 
energy loss process then determines the electron behavior depends mainly on the magnetic 
field strength and the size of the source. A simple comparison of the above relations shows 
that when Pad > -Psyn, i-e., 

Bitrmt) < n i = 2.3 x 10^^ n 1 (12) 

ar Vc/7 Vc/7 

adiabatic losses dominate over radiative losses. For example, if one considers expansion at 
speed u ^ c and an initial source dimension Rq ~ 10^^ cm, then for energies below 7 ~ 10'' 
the magnetic field can be as large as -B ~ 0.1 G and for energies less than 7 ~ 10^ the 
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adiabatic losses are still dominant for a value of 5 ~ 1 G. If the expansion of the source 
would not affect the hard slope at TeV energies, this could thus allow for a relaxation of the 
values used for SSC modeling of the source. In order to investigate this scenario in more 
detail, one needs to solve the electrons' kinetic equation 

?^^^(PM))-^^QM, (13) 

at (77 Te 

where Tg is the characteristic escape time and A^e the differential electron number. For 
simplicity, we neglect the escape term (rg oo), assuming that the sources expands with 
relativistic speeds u ~ 0.1c. For a constant expansion rate and continuous injection with 
rate Q(7, t) — )■ (5(7, R), we can replace the time variable t by the source dimension R . Then, 



the general solution of the kinetic equation (eq. (30) in lAtoyan fc Aharonianlll999l ) for the 
case of dominance of adiabatic losses is reduced to 

Neil, R) = ^No (^1) + - T-Q f-7, r) dr , (14) 
Ro \Ro J uj^ r \r J 

where the first term corresponds to the initial conditions, the contributions of which quickly 
disappears, and the second term relates to the continuous injection of relativistic electrons. 
Rq is the source dimension at the initial time to- 

We consider zero initial conditions [Nq = 0) and power-law injection of relativistic 
particles at constant rate 

Q(7, R) = Qol-"' 6(7 - 7o,min)e(7o,max " l^R - Ro) , (15) 
where 9 denotes the unit step function 

e(x-Xo) = l ^' ^^""^ (16) 

Pi > is the momentum index and Rq the radius at which injection starts. At radius R, 
electrons with initial cut-off energies 7o,max and 7o,max will have energies 7R,min and 7i?,,max, 
respectively, as they evolve according to Eq. ffTTl) . i.e. we have 

Moreover, there exists a critical radius i?* = -Ro7o,max/7o,min at which 7ij,max, i-e. the energy 
of the electron with initial injected energy 7o,max at R, becomes less than the initial 7o,min, 
so that the following two cases can be distinguished: 
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For i? < i?^,, or equivalently as long as 7R,max > 7o,min, we have 



iV(7, R) = { 7"^^ [1 - (f , 70,min < 7 < 7R,> 

PlU 



r 









Pi 


1 - 


( ' )"] 






\ 70, max / 



7R,max ^ 7 ^ 70,max 



-PI _ / Rr, 
'O.min I _Ro 



-Pi 



7R,min < 7 < 70,mir 



For R > R^:, or equivalently as long as 7R,max < 7o,min, the solution is 



Piu 



1 



-pi 



I O.min 



-PI _ / _R7 
/0,min I ifo 



-( ' 


)"] 


\ 70, max 




— Pi 

To, max' 


7r, 



, 70,mm < 7 < 70,n 



-Pl 



7R,min < 7 < 7R,max 



(19) 



The two solutions exhibit the same behavior. The differential electron number density drops 
with radius as ne{j,R) = '^oiaml ^ ^"^^ above the initial low-energy cut-o ff 7o,min 

adiabatic losses do not modify the power-law index ne(7, t) cx 7"^^ ( lKardashevlll962l ) . Below 
7o,min the resulting distribution is constant with respect to the electron energies, ne(7, R) cx 
7° (Fig. H]). The 7°-part of the electron population does not show up in the spectrum as 
the contribution of the injected 7o,min electrons (generating a 1/3- synchrotron wing) remains 
dominant at low energies (Figs. E] and E])- Thus, in contrast to the simple synchrotron 
cooling case, one has oc z/^/^ below the injected cut-off. For this reason, the classical hard 
spectrum picture at the TeV range can remain for timescales analogous to the source size. 
Even though electrons cool adiabatically as the source expands, the hard 1/3-synchrotron 
slope always appears below the synchrotron frequency related to the initial minimum Lorentz 
factor 

2 o.o^ (20) 



.syn 

min 



,BiR) cx — . 

1 \ ' pn 



Note that any decrease of this break energy occurs due to a decrease of the magnetic field. 
This is different to the pure synchrotron cooling case, where the corresponding break en- 
ergy follows the evolution of the minimum electron energy so that z/^" cx The same 
consideration holds for the energy regime where Compton scattering occurs. When we are 
deep in the KN regime (&(-R)7omin > 1) the energy below which the hard slope remains is, 
as mentioned above, 

(21) 



V, 



c 



« 7o,i 
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It therefore does not move to lower energies though the corresponding synchrotron frequency 
does. (In the pure synchrotron cooling case, obviously, z/^j^ oc 1/t^). As the source expands 
and the magnetic field drops, there will be an instant t corresponding to a radius R at which 
the KN regime no longer applies, and the break Compton frequency becomes 

^Zin « i^(^)7o,min ^ ^ . (22) 

which now moves to lower frequencies with the same rate as the synchrotron one. (Note that 
this reveals a very different time- dependence compared to the pure synchrotron cooling case 
where now z/^j^ oc 1/t^). However, the peak of the Compton flux still remains close to the 
initial 7o,min energy (Fig. [7]), and in total the decrease of the synchrotron peak flux is much 
stronger than the decrease in the Compton peak flux (Fig. [8]). In general, the dependence 
of the magnetic field on the radius R has important consequences for the behavior of the 
system even though the synchrotron losses are not important. The synchrotron peak flux 
varies as 

F^'^" oc NeB{Rf (23) 

and as we know from the solution of the kinetic equation that oc i?, the variability of the 
synchrotron luminosity should reflect the magnetic field dependence. The Compton flux, on 
the other hand, does not necessarily vary quadratically with respect to the synchrotron flux. 
As discussed above, the drop of the minimum Compton energy (which occurs naturally within 
the expansion-scenario) reduces the suppression of the cross-section and thereby supports 
the Compton emission. The variability pattern after "saturation" can therefore approach 
a quasi-linear dependence (cf. Fig [8]). Initially, during the raising phase before the two 
luminosities reach their maximum, the Compton flux can vary much more strongly, almost 
more than quadratically, with respect to the synchrotron flux. Moreover, close to saturation 
the Compton luminosity can exhibit a delay with respect to the synchrotron one as it reaches 
its maximum at later times compared to the synchrotron luminosity. 

The above considerations apply to situations where the expansion of the source com- 
pletely determines the evolution of the system. In reality, synchrotron losses could modify 
the electron distribution at high energies, namely for 

^ > ^- = (c) ■ '2^) 

However, as synchrotron losses decrease faster than adiabatic losses, one only needs to ensure 
that initially 7* > 7o,min- The change of the electron power index from —pi to — pi — 1 due to 
synchrotron cooling (cooling break) above 7o,min would then not disturb the hard 1/3 slope 
in the TeV range. 
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4. The external Compton case 



An alternative hypothesis to the SSC scenario concerns the Comptonization of a ra- 
diation field external to the electron source. In general, the optical-UV radiation field 
produced by a standard accretion disk could represent a non-negligible external source of 
photons to be up-scattered to the VHE 7-ray part of the spectrum. This radiation field 
could be up-scattered either directly by the rel ativistic electrons o f the jet (with target 
photons coming directly from the accretion disk, iDermer et al.l 119931 ) or more effectively af- 
ter being reprocess ed/re-scattered by emission line clouds like the broad line region (BLR) 
( ISikora et al.lll994j ). In external Compton (EC) scenarios, the geometry of the source and 
the location of the photon field with respect to the jet are of high importance as they can 
result in strong boosting or de-boosting effects on the photon energies. Here we explore the 
possibility of producing a hard TeV spectrum within the EC approach. We consider the BLR 
case, where the photon field is strongly boosted in the frame of the jet and up-scattered to 
higher energies. 

Let us consider a blob of relativistic electrons that travels with the jet of bulk Lorentz 
factor r along the z-axis. The jet passes through a region assumed to be filled with isotropic 
and homogeneous photons that obey a Planckian distribution of temperature T ~ 10^ — lO^i^ 
(corresponding peak frequency z/^ = 2.82kT/ /i ~ 5 x 10^^ — 10^^ Hz) . The central disk photon 
field is then characterized by a special intensity 



c [ekT 



-1) 



(25) 



a frac tion £ < 1 of whic h we assume is isotropized by re-scattering or reprocessing in the 
BLR ( ISikora et al.ll2002l ). so that the spectral energy density of the target photon field is 



BLR 



(26) 



where = An r^I^ is the spectral luminosity of the disk and where we take ~ for the 
disk radius. 

In order to take anisotropic effects into account, we transform the electron distribution 
from the comoving blob frame K' to the rest f rame K of the external phot on field, which 
in our case coincides with the observer's frame (IGeorganopoulos et al.l 1200 ll ). Electrons are 
assumed to be isotropic in the blob frame K'. In the photon frame K they exhibit a strong 
dependence on the angle 6, which is the observer's angle. As the up-scattered photons 
travel inside a cone I/7, we can make the approximation that they follow the direction 
of the electrons. The angle between the electron momentum and the bulk velocity of the 
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jet coincides with the observer's angle. The observer practically sees radiation only from 
electrons that in the photon frame are directed towards him. The observed flux then is 

F-. = ^ j Kij) W^(^,eph,e^) nph(eph)rf^ de^^ (27) 
where Nf,{E) denotes the differential number of electrons per energy per solid angle, and 



W = ,frL is the scattered photon spectrum per electron (given in iBlumenthal fc Gould 



19701 ). The unprimed quantities refer to the external photon field frame with number density 



^ph(eph) and the primed ones to the blob rest frame. 

We show the calculated Compton spectrum for a Maxwellian electron distribution in 
Fig. [9l The resulting TeV slope appears even harder than in the SSC case, with a limiting 
value of Fi, oc . Any photon field which is softer (flatter) than Fi, oc will dominate 
the Compton spectrum at low energies, as in our SSC model case where the up-scattered 
(synchrotron) photon spectrum follows F^ oc z/^^^. In all other cases, like in the external 
Compton scenario with a Planckian photon field (that at low energies follows F^ cx z/^), the 
characteristic behavior of the Compton cross-section appears, implying that the Compton 
spectrum at low energies (i.e., below ~ 7^60 where 7c is the electron break frequency in 
the photon rest frame) is dominated by the contribution from the up-scattering of the peak 
photons with ec ~ SfcT, yielding a. F^ oc u dependence. 

A similar consideration holds for a narrow (energetic) power-law electron distribution 
in an expanding source scenario. A hard VHE component F^ (x v should then appear below 
57o,min- The critical energy below which one can see this hard behavior of the Compton flux 
will not move to lower energies as the external target photon field is quasi-stable, so that 
the condition for the deep Klein-Nishina regime 7o,minec > 1 does not change. 



5. Summary 

The observed hard 7-ray spectra of TeV Blazars are difficult to explain within the most 
popular leptonic synchrotron-Compton models. The n{E) oc E'"^ shape, that the electron 
energy distribution is expected to quickly develop due to synchrotron cooling, usually results 
in a. Fy (X v^^l"^ radiation spectrum and therefore represents the limiting value of how hard 
the up-scattered spectrum can be in the TeV range. Moreover, modification due to Klein- 
Nishina effects can make the up-scattered TeV spectrum steeper than this and shift the 
Compton peak to much lower energies than observed. 

However, intrinsic source spectra as hard as Fy oc z/^^/^ exist, even when one only 
corrects for the lowest level of the EBL, most notably in the case of lES 1101-232 and 
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lES 0229+200 (lAharonian et al.ll2006l . l2007bl ). Most likely, the real source spectra are even 



harder. Investigating the possibility of forming hard VHE blazar spectra appears therefore 
particularly important. Methodologically, it seems necessary to first examine the "conven- 
tional" radiation and acceleration mechanisms, that have often been successful in interpreting 
Blazar observations, before adopting very different and often more extreme solutions. 

The results of this work show that within a simple homogeneous one-zone SSC approach, 
a power-law particle distribution with a large low-energy cutoff can in principle produce a 
hard {a = 1/3) - slope in the VHE domain {F^, oc z/") by reflecti ng the characteristic low - 



energy slope of the single particle synchrotron spectrum (cf. also iKatarzyhski et al.l 120071 ). 
As shown in section |U even harder VHE spectra approaching oc u [a = 1) can be achieved 
in the external Compton case for a Planckian-type ambient photon field. 

A power-law elec tron distribution with a high low-energy cut-off has been used in 



Tavecchio et al.l ( l2009l ) to model the emission from lES 0229+200 within a stationary SSC 
approach. In order to avoid the above noted synchrotron cooling problem, an unusually 
low value for the magnetic field strength was employed, leaving the particle distribution 
essentially unchanged on the timescales of several years. This goes along with a strong devi- 
ation from simple equipartition by several orders of magnitude (i.e., us/ue ~ 10~^). While 
it is known from detailed spectral and temporal SSC studies of the prominent 7-ray Blazar 
Mkn 501 that TeV sour ces may be out of equipartition at the one per- mill level or less 



( iKrawczynski et al.ll2002l ). the SSC modeling of lES 0229+200 suggest the hard spectrum 
sources to belong to the more extreme end. (Within external Compton models values closer 
to equipartition may be achieved, depending on the external photon field energy density). 
On the other hand, a large electron energy density (strongly exceeding the magnetic field 
one) could well facilitate an expansion of the source, and this motivates a time-dependent 
analysis: 

Using a time-dependent SSC model, we have shown that the hard (a = 1/3)- VHE slope 
can be recovered, when adiabatic losses dominate over the synchrotron losses for the low- 
energy part of the electron distribution (i.e., for Lorentz factors less than the injected 7min)- 
The main reason for this is, that the resultant electron distribution below 7min becomes flat 
and therefore does not show up in the SSC spectrum. Interestingly, this scheme also allows 
one to relax the very low magnetic field constraints. 

We also examined the relevance of a Maxwellian-like electron distribution that peaks at 
high electron Lorentz factors ~ 10^. Such a distribution represents a simple time-dependent 
solution that already takes radiative energy losses into account, and turns out to be capable 
of successfully reproducing the hard spectra in the TeV range (with limiting values a = 1/3 
and a = 1, respectively). Maxwellian distributions can be the outcome of a stochastic 



- 16 - 



acceleration process balanced by synchrotron or Thomson cooling. Depending on the physical 
conditions within a source, e.g., if particles undergo additional cooling in an area different 
from the acceleration one (jSauge &: Henri! l2006l : iGiebels et al.l 120071 ) , or if the medium is 
clumpy supporting a "multi-blob" scenario in which the observed radiation is the result of 
superposition of regions characterized by different parameters, the combination of pile-up 
distributions may allow a suitable interpretation of different type of sources. For the case 
presented here, they demonstrate a physical way of achieving the high low-energy cut-offs 
needed in leptonic synchrotron-Compton models for the hard spectrum sources. 

Although our main purpose here is not to fit data. Fig. [10] shows that a Maxwellian- 
type electron distribution could also provide a satisfactory explanation for the hard TeV 
component in lES 0229+200. 

Our results illustrate that even within a leptonic synchrotron-Compton approach rela- 
tively hard intrinsic TeV source spectra may be encountered under a variety of conditions. 
While this may be reassuring, the possibility of having such hard source spectra within 
"standard models" unfortunately constrain the potential of extracting limits on the EBL 
density based on 7-ray observations of Blazars, one of the hot topics currently discussed in 
the context of next generation VHE instruments. 
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Fig. 1. — Stationary SSC spectra for different values of the low-energy cut-off. Above 
7min ~ 3 X 10^ we are very deep in the Klein-Nishina regime and the peak of the Compton 
emission appears very sharp. As one reduces 7min5 the suppression of the cross-section 
decreases and the minimum Compton energy drops to lower energies. Thus, the peak of 
the Compton flux raises significantly, whereas the synchrotron peak remains constant. A 
Doppler factor 6 = 50 has been used for the plot. 
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Fig. 2. — Evolution of the observed SSC spectrum for constant injection of a narrow power- 
law electron distribution, with modifications due to synchrotron cooling taken into account. 
The magnetic field is B — IG. The hard (1/3) synchrotron and Compton spectral wings 
are observed for timescales shorter than the cooling timescale of the 70-particles, i.e., in 
the present application for timescales < 0.1 days. The figure shows the expected spectral 
evolution for a total (observed) time t ~ 1 day. Parameters used are Rq = 7.5 x 10^^ cm. 



7n 



7 X 10^ 7„ 



2 X 10^, power law index p — 2.85 and Doppler factor S — 25. The 



total injected power is Q ~ 10 erg/sec. 
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Fig. 3. — SSC modelling with different electron distributions. Black line: Power-law 
with large value of the minimum energy (as in Tavecchio et al. 2009). The parameters 
used are 7mm = 5 x 10^, 7max = 4 x 10^, power law index p = 2.85, B = A x lO"'' G, 
fee = 6.7 X 10* cm~^, R = 5.4 x 10^^ cm and Doppler factor 5 = 50. Red line: Relativistic 
Maxwellian distribution N^. = Kej"^ exp{ — — ) with parameters 7c = 1.5 x 10^, B = 0.07 G, 
Ke = 3 X 10~^^cm~^, R = 2 X 10^'^cm and 6 = 33. The peak of Compton flux occurs in 
the KN regime as {B/Bcr)^^ — ISO >> 1. Blue line: Relativistic Maxwellian distribution 
Ne = Ke-f^exp^-^J^ with parameters 7^ = 5.3 x 10^ B = 0.06 G, = 4 x lO'^^ cm'^, 
R = 2x 10^^ cm and 6 = 33. 
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Fig. 4. — Illustration of the evolution of the electron distribution for constant injection 
and dominant adiabatic losses. The expansion of the source does not modify the power-law 
index above the initial low-energy cut-off 7o,min, whereas below it the distribution becomes 
approximately flat. The electron number density depends on radius as Ue oc 
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Fig. 5. — Evolution of the observed SSC spectrum with constant injection of a narrow power- 
law and dominant adiabatic losses. As the source evolves, the synchrotron peak decreases and 
gets shifted to smaller energies following the decrease of the magnetic field. For the magnetic 
field, we use an initial value Bq = 0.075 G and we assume that it scales as B = Bq (Rq/R) 
(i.e., m = 1). The initial radius is Rq = 7.5 x 10^^ cm, expanding up to i? = 10i?o (at 
u = 0.1 c) and corresponding to observed timescales of the order of t ~ 30/5 days. The 
total injected power is Q ~ 5 x 10^^ erg/sec. Other parameters used are 7min = 3 x 10^, 
7max = 2 X 10'', power law index pi = 2.85, Doppler factor 6 = 25 and Qo = 1.5 x 10^^ sec^^. 
Note that timescales are comparable to the synchrotron cooling case. 




Fig. 6. — Same as Fig. [5] but for a different magnetic field scaling, m = 2. Other parameters 
are the same as in Fig. O 
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Fig. 7. — Evolution of the SSC spectrum with dominant adiabatic losses (for B oc 1/R, i.e., 
m = 1) from Rq = 7.5 x 10^^ cm to R = 10i?0; and corresponding observed luminosities. 
Whereas the synchrotron peak gets with time significantly shifted to lower energies, the 
Compton peak can appear almost static. The Compton flux reaches its maximum at a 
greater radius than for the synchrotron one. Synchrotron and Compton fluxes are shown 
until maximum (with red lines) and after (with black lines). 
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Fig. 8. — Evolution of the SSC spectrum luminosities with dominant adiabatic losses (for 
m = 1), see also Fig. [71 The Compton flux reaches its maximum at greater radius (i.e., 
later) compared to the synchrotron one. While during the raising phase the variability 
pattern approximately shows a quadratic behavior, the correlation becomes almost linear 
during the declining phase. 
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Fig. 9. — External Compton scenario for a Maxwellian-type electron distribution (with 
Op = 0). The observed spectrum is calculated for different angles 9 to the observer. The 
synchrotron slope follows oc u^/^ . In the TeV range (x ^ i.e., harder than in the SSC 
case. The dashed line corresponds to the assumed disk spectrum. The bulk Lorentz factor 
of the jet is r = 13 and the peak energy of the electron distribution is 7c = 2 x 10^. For 
the disk photon field a temperature T = 1.75 x 10"^ K is assumed. The relevant radius Rd of 
the disk is considered to be of the same dimension as the jet (10^^ cm). The magnetic field 
is 5 = 1 G and a fraction ^ = 0.1 of the disk photons is assumed to be rescattered by the 
BLR. 
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Fig. 10.— The hard spectrum blazar lES 0229+200 at z=0.139 with SED modeled within 
an SSC approach using Maxwelhan-type electron distributior is. All parameters u s ed ar e 
the same as in Figj3l Data points shown in the figure are from IZacharopoulou et al.l ( 1201 ll ). 
w here the intrinsic (d e -abso rbed) source spectrum has been derived based on the EBL model 
of iFranceschini et al.l (120081 ) with (i) EBL level as in their original paper ("low level EBL") 
and (ii) (maximum) EBL level scaled up by a factor of 1.6 ("high level EBL"). 



